New Parametrization for the Scale Dependent Growth Function in General Relativity 
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We study the scale dependent evolution of the growth function 8(a, k) of cosmological perturba- 
tions in dark energy models based on General Relativity. This scale dependence is more prominent 
on cosmological scales of lOO/i" 1 Mpc or larger. We derive a new scale dependent parametriza- 
tion which generalizes the well known Newtonian approximation result fo(a) = d }? g ° = fi(a) 7 
(7 = TT f° r ACDM ) which is a good approximation on scales less than 50/i _1 Mpc. Our gener- 

-M^- where f(o,fc) = ^S 22 



alized parametrization is of the form f(a) 



i+S(a,k) s^~;- ak 2 ■ We demonstrate 

that this parametrization fits the exact result of a full general relativistic evaluation of the growth 
function up to horizon scales for both ACDM and dynamical dark energy. In contrast, the scale 
independent parametrization does not provide a good fit on scales beyond 5% of the horizon scale 
(k ~ OMh^Mpc). 



PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 
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INTRODUCTION 



The observable growth function 5(a) = -^(o) of the 
linear matter density contrast as a function of the cos- 
mic scale factor a (or equivalently the redshift z) provides 
a useful tool to test theoretical models attempting to ex- 
plain the observed accelerating expansion of the universe. 
In particular, the combination of the observed expansion 
rate of the universe H (a) [3, [H, Q with the growth func- 
tion 6(a) 0, i, i, 0, 1, 1 can provide significant insight 
into the properties of dark energy [lCJ, LLI[ driving the 
accelerating expansion (e.g. sound speed, existence of 
anisotropic stress etc) or even distinguish it from modi- 
fied gravity [l2| which may alternatively be responsible 
for the accelerating expansion. 

The sub-Hubble (scale independent) growth function 
can be obtained for any homogeneous dark energy cos- 
mology in the context of General Relativity by solving 
numerically the growth equation On larger scales, 

there are several numerical studies solving the full cou- 
pled relativistic equations for cosmological perturbations 
and making very precise predictions about structure at 
large-scales and/or high redshift (eg 14]). Alternatively, 
the growth of cosmological perturbations may be approx- 
imated by using a simple parametrization. The standard 
parametrization of the linear growth function 6(a) is usu- 
ally made by introducing a growth index 7 defined by 



fo(a) 



din Sq 



Sl m (a) n 



where a 



d\na 

is the scale factor and 



(1.1) 



(1.2) 



is the ratio of the matter density to the critical density 
when universe has scale-factor a where Hq is Hubble con- 
stant and flam is ratio of mass density to critical den- 
sity. This parametrization [l5| provides an excellent fit 



to the evolution equation for 6(a) in the small scale (sub- 
Hubble) approximation 



2H6 - AnGp m 5 = 



(1.3) 



where an overdot denotes the derivative with respect to 
time and p m is the matter density. Changing variables 
from t to In a we obtain the evolution equation for the 
growth factor / as 



H 



f' + f + fi-jp+Z) 



(1.4) 



where ' = d/dlna. For dark energy models in a flat 
universe with a slowly varying equation of state w(a) = 
= Wo, the solution of eq. (|1.4[) is well approximated 
by eq. (jO) with [3 



3K - 1) 
6wq — 5 



(1.5) 



which reduces to 7 = yj for the ACDM case (wq = — 1). 
It is therefore clear that the observational determination 
of the growth index 7 can be used to test ACDM [4] . It 
has been shown Q that even in the context of dynamical 
dark energy models consistent with Type la supernovae 
(Snla) observations the parameter 7 does not vary by 
more than 5% from its ACDM value. However, in the 
context of modified gravity models 7 can vary by as much 
as 30% (e.g. for the DGP model^ j DGP ~ 0.68 0) 
while scale dependence is also usually introduced @, [13, 

El. 

Current observational constraints on 7 are based on 
redshift distortions of galaxy power spectra [17| , the rms 
mass fluctuation <j%(z) inferred from galaxy and Ly — a 
surveys at various redshifts [H!]-[l||, weak lensing statis- 
tics [10] , baryon acoustic oscillations [2l| , X-ray luminous 
galaxy clusters [22| , Integrated Sachs- Wolfe (ISW) effect 
|23| | etc. Unfortunately, the currently available data are 
limited in number and accuracy and come mainly from 
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the first two categories. They involve significant error 
bars and non-trivial assumptions that hinder a reliable 
determination of 7. Thus, the current constraints on 7 
are fairly weak [1] and are expressed as 



v = 674+ - 195 

7 — u -°'*_0.169 



(1.6) 



This however is expected to change in the next few years 
when more detailed weak lensing surveys are anticipated 
to narrow significantly the above range. 

A crucial assumption made in the derivation of eq. 
(|1.3|) in the context of general relativity metric pertur- 
bations is the assumption that the scale of the perturba- 
tions is significantly smaller than the Hubble scale 
This assumption however does not lead to a good ap- 
proximation on scales larger than about 50h~ 1 Mpc [241 ] . 
In fact, the perturbed metric of spacetime takes the form 
(in the Newtonian gauge): 



ds 2 = -(l + 2$)d< 2 



(1 - 2<$>)a 2 lij dx l dx 3 , 



(1.7) 



where 7^ is the metric of the spatial section and we are 
ignoring anisotropic stresses. The evolution of density 
perturbations on all scales is dictated by combining the 
background equations 



H 2 
P 



8nG 



(Pm + Pde) 



-3H(p + p) 



(1.8) 
(1.9) 



(assuming a flat universe with only pressureless dark 
matter and (non-clustering) dark energy) with the per- 
turbed linear order Einstein equations in the Newtonian 
gauge 25] (p m and p^e & re the matter and dark en- 
ergy densities respectively while p = wp is the pressure). 
The resulting (anisotropic stress-free) equations are of 
the form 



$ = -4ff$ + 8irGp de w de $ 
k 2 

5 = 3$ + —v f 
a 2 



with constraint equations 



3H (H$ + $) 



(F$ + <E>) 



-AnGSpm 
-4irGp m Vf 



(1.10) 
(1.11) 
(1.12) 

(1.13) 
(1.14) 



where $ is the Newtonian potential, Vf = —va (v is the 
velocity potential for dark matter) and we have general- 
ized the derivation of Ref. [24| to the case of a general 
dark energy equation of state parameter Wd e (a)- Clearly, 
equations (|1.10[) - (|1.12p involve a scale k dependence in 
contrast to the small scale approximate equation (|1.3p 
which is scale independent. It has been demonstrated 
[IH that for Wde = -1 (ACDM ), the solution of equa- 
tions (|1.10p - (|1.12p (with initial conditions for $ and vt 
derived from the constraint equations with <f> = 0) devi- 
ates significantly from the solution of the approximate 



1.0 Scale independent parametriztitigiL——^^ 



Scale Dependent G.R. Solution 




Scale dependent parametrization 



JO 









Scale independent parametrization 


/s~ Scale Dependent G.R. Solution 


k=0.0O4h Mpc~' 


// 




// 
■ 

1 




I 

! Scale dependent parametrization 

/ 

/ 




2 4 


6 8 10 



Scale independent parametrization 



Scale Dependent G.R. Solution k=ojm t> Mpc-' 
/■/ Scale dependent parametrization * v, "*>» 



FIG. 1: a: The growth factor / obtained from the solution of 
the general relativistic system (k = O.Ol/i -1 Mpc, Qom = 0.3, 
ACDM , dotted line) compared with the scale independent 
parametrization (continuous line) and the corresponding gen- 
eralized scale dependent parametrization (thick dashed line) 
b: Similar to a. for k — 0.004/i -1 Mpc. c: Similar to a. for 
k = 0.001/i _1 A'/pc. 



equation (|1.3p on scales larger than about 50h~ 1 Mpc 
(k < 0.02hMpc). This is demonstrated in Fig. 1 where 
we compare the growth factor f(z) based on the solu- 
tion of the full general relativistic linearized equations 
(|1.10P - (|1.12p with the corresponding growth factor of the 
standard parametrization (jl.ip for k = O.OOlhMpc -1 , 
k = 0.004/iMpc" 1 and k = O.OlhMpc- 1 (w de = -1 and 



fiom = 0.3). 
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FIG. 2: a: The accuracy of the scale independent 
parametrization in terms of A as defined in Eq. (|2.16|l . b: 
Likewise for the scale dependent parameterization. 



Since this range of scales is accessible to present and 
especially to future large scale structure observations, 
it becomes clear that there is a need for an improved 
parametrization for the growth factor / that will be com- 
patible with the full general relativistic solution up to 
horizon scales. The derivation of such a parametrization 
is the goal of the present study 1 . 

The structure of this paper is as follows: In the 
next section we focus on the case of ACDM and show 
that there is a scale dependent generalization of the 
parametrization (jl.ip that approximates well the full gen- 
eral relativistic linear growth solution up to horizon scales 
(k ~ 10 hMpc ). In section III we generalize this re- 
sult to the case of dynamical dark energy models and 
finally in section IV we conclude and point out future 
extensions of this work. 



1 Note that if one works in the synchronous gauge, then Eq. dl .31 
is exact (see e.g. the discussion in [241]'). However, the focus of 
this work is the Newtonian gauge. 



2. SCALE DEPENDENT GROWTH 
PARAMETRIZATION IN ACDM 

We start by sketching the derivation of equation (jl.3j) 
in the context of general relativity using the sub-Hubble 
approximation. The linear matter overdensity 5p m may 
be expressed p5j in terms of the gravitational potential 
<& and the background variables as follows: 



k 2 

4nG5p = 



3i? 2 $ 



(2.1) 



In the sub-Hubble (small scale) approximation (^ 3> 
H 2 ) equation (|2.ip takes the form 



k 2 

inGSp = — <E> 
a 1 



(2.2) 



where we have also assumed a slowly varying gravita- 
tional potential On dimensional grounds we would 
expect that 37J 2 <I> ~ 3i?$ and thus it may seem un- 
natural to drop the 3if<I> while keeping the 3_ff 2< f>. For 
example the effective combination of both terms would 
be 6H 2 &. However, as discussed in section 3 where an 
arbitrary coefficient is considered for the term H 2 Q, the 
general relativistic solution is best approximated by a 
term of the form 3H 2 Q which implies that the time vari- 
ation of the gravitational potential $ can be ignored. 

The general relativistic equations (|1 . 10|) - (| L . 14|) lead to 
the following equation for the matter overdensity S 



k 2 

6 + 2H6+—$ = Q 



(2.3) 



Using the sub-Hubble approximation (|2.2p we obtain the 
scale independent approximate equation (|1.3|) . On the 
other hand, if we avoid this approximation in equation 
(|2.1[) . solve for $ and substitute in equation (|2.3[) we 
obtain the following scale dependent evolution equation 
fortfO: 



where 



(2.4) 



(2.5) 



The solution of equation (|2.4p provides a much better ap- 
proximation to the full linear general relativistic system 
(|1.10[) - (|1.12[) up to horizon scales. On scales larger than 
the horizon even equation (|2 .4[) breaks down since on 
these scales, the time derivative of $ can not be ignored. 

Given the successful approximation of the solution of 
(|2.4p to the exact linear general relativistic solution, 
it becomes important to construct a scale dependent 
parametrization that is analogous to (jl.ip and solves (ap- 
proximately) (|2.4p for all scales k. In order to construct 
such a parametrization we focus on the matter dominated 
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FIG. 3: a: The values of the parameters (7, a) that fit best the general relativistic growth solution as functions of the dark 
energy equation of state parameter wq (wi was fixed to 0). Blue disks correspond to a, black triangles to 7. The error bars in 
a describe the la variation of a as the scale k is varied in the range k g [0.001, O.ljft" 1 Mpc b: Similar to a for the dark energy 
equation of state parameter wi (wo was fixed to -1.4). 



era when most of the growth occurs and express £(a, k) 
as 



£(o,k) = 



111 

ak 2 



(2.6) 



Equation (|2.4p may be expressed in terms of the growth 
factor f — in the form 

J a In a 

/' + f + Ho„(«,)HTT!i& <"> 

where ' = . and we have assumed ACDM for Hid). 

a In a ' V ' 

Notice that the matter domination approximation is im- 
plemented explicitly only with respect to the scale de- 
pendent correction. The sub-Hubble scale independent 
parametrization does not explicitly set fi m (a) = 1 but 
its derivation implicitly assumes £l m la) ~ 1 (see eg 
[H[). Thus we have used f2 m (a) in eq. (|2.7|) in order 
to make connection with the sub-Hubble scale indepen- 
dent parametrization fo(d). 

For sub-Hubble scales £(fc, a) — > and equation (|2.7p 
reduces to (| 1 .4[) whose solution is well approximated by 
(PnH with 7 = p-. 

We now consider a small £ and look perturbatively for 
a solution of the form 

/(M) = / (a)(l-/i(M)) (2-8) 

Substituting (|2.8p in (|2.7[) and keeping only terms linear 
in fi and £, it is straightforward to show that 

111 



/i(fc,a) = £{k,a) = 



ak 2 



(2.9) 



where we have also used the fact that in a matter domi- 
nated universe 5(a) ~ a. 

This analysis may be extended to arbitrary order in 
perturbation theory by expanding (1 + £(fc, a)) -1 



N 



(l + e(fc,a))- 1 -^(-l)"e(fc,a) n 



(2.10) 



and looking for a solution of the form 

N 

/(fc,a) = /o(a)^(-l)"/«(fc,a) r 



(2.11) 



It is then straightforward to show order by order that 
f n (k,a) = £(fc,a). Therefore, the parametrization 



f{k,d) 



Ma) 



fl m (a)"' 



1 + Z(k,d) 1 + 



(2.12) 



is an approximate solution of equation (|2.7[) and provides 
a good approximation to the solution of the general rel- 
ativistic system (|1.10p - (|1.12p up to horizon scales. 

Alternatively, the form of f(k,a) for a ACDM back- 
ground can also be deduced in a manner similar to the 
appendix in [l5| . Guided by the matter-dominated case, 
we assume that 



f(k,a) 



ak 2 



1 



(2.13) 



where X(a) is some arbitrary time-dependent function. 
We next switch to x = 1 — il m (a) as the time variable, 
assume the ACDM evolution for f2 m (a), and express the 
left hand side of Eq. (|2.7p as a series in x: 



,1/3 



k 2 X - 3H 2 k 2 n m0 
2Hl (1 - il m0 ) 1/3 ftfloX, 



k*X 2 - 9H^Q 2 n0 - 15H^Pnl X 1 \ 2 



6H*(i-n m0 ) 2 / 3 n^x 2 



x 2 / 3 + O [x] 



(2.14) 



where we have used a Taylor expansion for the function 

X(x) 



X(x) = Xo+Xtx + O [x 2 ] 



(2.15) 
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For Eq. (|2.13|) to be an approximate solution to 
Eq. (|2.7[) . we can expect the leading order coefficients 
of the expansion (|2. 14|) to vanish. The vanishing of 
the first two coefficients immediately implies that that 
Xq = 3iJof2 TO o and X\ = 0. The latter fact makes it 
clear that any x dependence of X is weak (quadratic or 
higher). (Note that higher coefficients of (|2.14|) contain 
non-linear combinations of the coefficients of X(x), mak- 
ing it difficult to continue this process to deduce further 
coefficients in Eq. (j2.15p .) We therefore conclude that 
Eq. (|2.12[) provides a good approximation to the solution 

of Eq. am . 

The accuracy of our parameterization is demonstrated 
in Fig. 1 where we compare the form of the scale de- 
pendent parametrization Eq. (I2.12|) to the general rela- 
tivistic numerical solution with the corresponding fit of 
the standard parametrization for three different scales. 
It is clear that up to approximatelly the Hubble scale 
(k ~ O.OOlhMpc" 1 ) is accurate at a level better than 5% 
at least up a to redshift z — 10. We further quantify 
the accuracy of either parameterization through A, the 
fractional deviation of / predicted by the parameteriza- 
tion (i.e. Eq. for the scale independent case and 
Eq. (|2.12|) for the scale dependent case) from the exact 
/: 

^ /exact /parametrization ^ -j^g^ 

/exact 

In Fig. |a we display A for the scale independent 
parameterization for the same scales as considered in 
Fig. Jl}. This can be contrasted with the correspond- 
ing A for the scale-dependent case shown in Fig. [2]d. 

The analysis of this section has assumed a background 
expansion based on ACDM . In the next section we 
demonstrate that the parametrization (|2.12[) is also a 
good approximation to the general relativistic system so- 
lution in the case of a dynamical dark energy background. 



to the numerical solution and determine the best fit val- 
ues of the parameters (7, a) for various values of the scale 
k and of the dark energy parameters (wq,wi). In par- 
ticular, we first fix (wo,wi) and also the scale Ho a 
value much smaller than the horizon (eg k ~ lhMpc^ 1 ). 
We then solve numerically the system (|1.10|) - (|1.12| . find 
f(k 7 a) and fit the parametrization (|3.2p to the numeri- 
cal solution thus determining the best fit value of 7 (the 
value of a is irrelevant for large k). Next, we consider 
larger scales (up to k = lO^hMpc^ 1 ) and for each value 
of k we find the best fit value of a for the fixed value of 
7 obtained from the sub-Hubble fit. As anticipated from 
the discussion of the previous section we find that the 
best fit value of a is practically independent of the scale 
k and is equal to 1. Finally, we repeat the above steps for 
different values of (wq , w\ ) to find the dependence of the 
best fit parameters (7, a) on the dark energy properties. 

We find that for values of (wo, w±) consistent with ob- 
servations (7, a) ~ (0.55,1). This is demonstrated in 
Fig. 3 where we show the dependence of the best fit pa- 
rameter values for (7, a) in terms of Wg (for fixed Wi = 
in Fig. 3a) and in terms of w\ (for fixed wo — —1-4 in 
Fig. 3b). The error bars in the parameter a describe 
the la variation of a with the scale k. This is clearly 
negligible except of the pair (wo, w\) — (—1.4, 1.2) which 
implies significant presence of dark energy at early times 
thus affecting our assumption of matter domination at 
early times. 

In Fig. (H)) we plot A for the scale dependent and in- 
dependent parameterizations for the values (wq^wi) = 
(-1.4, 0) (top frame) and (w , Wi) = (-1.4, 1.2) (bottom 
frame) . 



3. SCALE DEPENDENT GROWTH 
PARAMETRIZATION WITH DYNAMICAL 
DARK ENERGY 

It is straightforward to solve the general relativistic 
system (|1.10j) - (|1.12p in an arbitrary homogeneous dark 
energy background i.e., in what follows, we ignore per- 
turbations of the dark energy component. However for 
a discussion of the growth of cosmological perturbations 
in which the dark energy component is not assumed to 
be homogeneous, see e.g. [26|, [27], [28| and references 
therein. For the evolution of iu, we consider a dark energy 
parametrization (29j 

w de {a) — Wq + Wi(l - a) (3.1) 

and solve the system (jl . 10|) - (II . 12[) for various values of 
the parameters (u>o, wi). We then fit the parametrization 



f(k,a) = 



Q m (a)i 



(3.2) 



CONCLUSION - DISCUSSION 



We have found a scale dependent parametrization of 



the growth function / 



dlni5 
d In a 



that is free from the sub- 



Hubble approximation of the standard parametrization 
(jl.ip . Our parametrization described by equation (|2.12p 
approximates very well the full linear general relativistic 
solution up to horizon scales and differs significantly from 
the standard parametrization on scales larger than about 
bOh^Mpc. 

An important extension of our work is the derivation of 
the best fit parameter values (7, a) in the context of mod- 
ified gravity models. It is also interesting to investigate 
whether a more general scale dependent parametrization 
is required to describe the growth function on large scales 
in modified gravity theories. 

The mathematica files used for the produc- 
tion of the figures may be downloaded from 



http: / /leandros.physics.uoi.gr/newpar.zip| 
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FIG. 4: A for the scale independent and scale dependent pa- 
rameterizations at a scale of k=.004h Mpc -1 for the (wo, wi) 
choices shown. 
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